EXPLICIT SOLUTIONS FOR THE EXIT PROBLEM FOR A 
CLASS OF LEVY PROCESSES. APPLICATIONS TO THE 
PRICING OF DOUBLE BARRIER OPTIONS 
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Abstract. Lewis and Mordecki have computed the Wiener-Hopf factoriza- 
tion of a Levy process whose restriction on ]0, +oo[ of their Levy measure 
has a rational Laplace transform. That allows to compute the distribution of 

(Xt, inf X s ). For the same class of Levy processes, we compute the distribu- 

0<s<t 

tion of (Xt , inf X s , sup X s ) and also the behavior of this triple at certain 

0<s<t 0<s<t 

stopping time, like the first exit time of an interval containing the origin. Some 
applications to the pricing of double barrier options with or without rebate are 
evocated. 



1. Introduction 

There are very few examples of Levy process for which the so-called " exit problem" 
can be explicitly solved (see [KK05] , |R90] . |R72j ). We present here this explicit 
solution for a class of Levy processes which has been introduced by Lewis A.L. 
and Mordecki E. [LM08 ]; that is the class of the Levy processes whose restriction on 
]0, +oo[ of their Levy measure has a rational Laplace transform. This happens when 
this restriction is a finite linear combination of exponential or gamma distributions. 
Lewis A.L and Mordecki E. [LM08] (see also Asmussen, S. Avram, F. and Pistorius 
M.R. [AAP04] ) have computed the so-called Wiener-Hopf factorization of these 
Levy process. That brings to closed forms for the distribution of the maximum of 
the process before an independent exponential time and for the joint distribution 
of (T x , Xtx) where T x is the first time where the Levy process X crosses upward 
a level x and Xt* its position at that time. 

As an application, when adopting the exponential Levy model, Y t — Y$.e Xt , for 
a financial asset Y t with X t of the preceding form (see for example |CT03j ). an 
immediate consequence of the preceding results is the computation of the (temporal 
Laplace transform) price of the double barrier option with this underlying asset. 
That is, an european option which is activated (in) or desactivated (out) when the 
asset {Yt)te[o,T] cross up (down) a barrier H before the time of maturity T. This 
application follows the computation of the price of the simple barrier option, by 
the Wiener-Hopf factorization (see |AAP04| ). 

In this paper, for the same class of Levy processes, we solve the exit problem. 
More precisely, we give closed form of the joint distribution of the minimum and 
maximum of the process before an independent exponential time and among other 
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behaviors of the process at certain stopping times, we give the Laplace transform 
of the joint distribution of (T^,X T b) where is the first time that the process 
leaves a bounded interval [—a, b] containing the origin and X T b is its position at 
that time. 

As an application, we mention how one can deduce the price of the double barrier 
option, in or out, that means that the option is activated or desactived if the asset 
crosses up a barrier b or down a barrier a in the exponential Levy model. 
The paper is organized as follow : In section 2, we recall the general results 
on Wiener-Hopf factorization, we give a shortened proof of the results of Lewis- 
Mordecki and settle few other preliminary results. For that, we use exclusively 
elementary complex analysis arguments. After introducing few more notations in 
section 3, we give in section 4, 5, 6, 7 and 8, all the results on the fluctuations of 
our class of Levy processes. In section 9, we recall the main result of |F10j which 
is the main tool of this work. The proofs follow in section 10, 11, 12, 13. 
These results, in the symmetric case, are related to what is called "Bargmann 
equations" in the litterature on inverse problems of the spectral theory (see for 
example |F63] V 

2. The assumption and the Wiener-Hopf factorization 

We suppose that A is a real Levy process possibly killed at an independent expo- 

nantial time and we denote £ the life time of A. 

Let <j> be the Levy exponent of A, so that the identity 

E( e - iuX *l t<c ) = e~*^ u ) 

is fulfilled for every time t and every imaginary number iu G zR, 4> is continuous 
on iR and 4>(0) is the rate of the exponential distribution of the life time £, 

c/)(0) = if and only if ( = +oo a.s. ("A does not die"). 

We now work under Lewis-Mordecki's assumption. It is based on the following 
rather obvious fact. 

Proposition 2.1. The conditions below are equivalent 

(i) The Levy measure of X, it, is of the following form 

n „ . 

n(dy) = y Cj — -e JiV dy on ]0,+oo[ 

J = l J 

n € N, rij E N, K(7j) > 0, Cj e C {j = 1, . . . , n) 

(ii) The exponent <f> is of the form 

P 

<j)(iu) = cj) (iu) + g^u), 

where <f)~ is the exponent of a Levy process without positive jumps, P and Q are 
polynomials and Q has all its roots on the complex half plane {A; 5R(A) < 0}. 

Proof. Assume (i). Notice that the compound Poisson process with Levy measure 

3=1 ny 
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has Levy exponent 

V ,_1 1 1 = PW 

The function <p — ^ is then the Levy exponent of a Levy process without positive 
jumps. This establishes (ii). 

Conservely, assume (ii). Since <p and <f>~ are exponents of Levy processes, we 
know that cp(iu) = 0(u 2 ) and <fi~(iu) — 0(u 2 ) (see proposition 2 chapter 1 of 
|B96| ). thus ^(iu), construed as a rational function of u, is 0(u 2 ). Consequently 
degP < degQ + 2. 

Now, write ^(iu) in its fractional expansion 

P n k- 

— (iu) = au 2 + bu + c + > — — — r 

Q 4^ (iU + 1j) n i +1 

Then ^ (iu) is the Fourier transform of the Schwartz distribution (here, 5 stands 
for the Dirac mass at 0); 



aS" - ibS' + cS + V k,—e-^ x l x>0 dx 

3 = 1 3 

On the other hand, one can deduce from the Levy Kinchin formula applied to cf> 
that <j)(iu) is the Fourier transform of a distribution whose restriction to ]0,+oo[ 
is — l x> QTr(dx) (see chapter 5 of |V02] for example) where 7r is the Levy measure. 
From this, we see that 



^x>o^(dx) = - ^2 — re llX l x> odx 
j=l nj " 

This finishes the proof. 

□ 

The following is a famous result (see proposition 2, chapter I of |B96j ) which applies 
to any Levy process and which will be of use later . 

Lemma 2.2. There exists an exponent of a subordinator tp and an exponent of the 
opposite of a subordinator ip, such that 

ip.Tp = (j) on iR 
Such a couple (if),ip) is unique up to a multiplicative constant. 

The functions ip and ip will be refered to as the "positive" (for ip) and "negative" 
(for ip) Wiener -Hopf factors of the exponent <fr, unlike the ordinary uses which 
affect these expressions to the functions and ^jp. 

ASSUMPTION : In the rest of this paper, conditions (i) and (ii) of 
proposition 12. 11 will be assumed to be satisfied and this assumption will 
be call Assumption 12.11 . 

Property (ii) of proposition 12.11 and the fact that (p~ has an holomorphic exten- 
sion on the half plane {5R(A) < 0} (see chapter 7 of |B96j ) imply that <j> has a 
meromorphic extension on this half plane, and we we will still denote by <f> this 
extension. 
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We will denote 

-71) -72, -In 

the poles of (p repeated according to their multiplicity. In particular, deg Q = n 
when the fraction ^ is irreductible, in the expression of condition (ii) of proposition 

12.11 {<p = <p~ + ^) and we will denote 

nj := ${k < j;j k =7j'} j = l,...,n 
Notice that the rij are all zero if and only if all the poles of <p> on the half plane 
{Re{\) < 0} are simple, or equivalently, if and only if the Levy measure is a (not 
necessarily positive) combination of exponential distributions. 

Theorem 2.3. ip is of the form : 

nr(A + A) 
^ (A) = ^n«(A + 7j )' 

where i/'oo is a positive constant. 

A 



//lim,\->-oo g]0,+oo] then m = n + 1. Otherwise, lim,\->.-oo G] — 00, 0] 



and m = n. 

More over, {— /3\, . . . , — f3 n } is the set of the roots of <j) lying on the half plane 
{5R(A) < 0} (repeated according to their multiplicity) together with if 4>(0) = 
and <£'(0) > 0. 

Proof. Notice first that, thanks to Levy Kintchin formula, 4>(X) does not vanish for 
A on the imaginary axis except possibly for A = 0, <f> is a meromorphic funtion on 
{5R(A) < 0} and 4> is continuous by the left at each point of the imaginary axis zR. 
On the other hand, the function %j> is an exponent of a subordinator, thus it is 
analytic on the open half plane {5£(A) > 0}and continuous on the closed half plane 
{5R(A) > 0} (see [B96] chapter 3 for example). Also, the identity ipip = <fi on zR 
implies that ip — on iR. Thus, since <ft is meromorphic and ip holomorphic on 
{5R(A) < 0}, we see that ip has a meromorphic extension on {5R(A) < 0}.We denote 
again ip this extension. 

More over ip = i is continuous by the left on the axis iR, except possibly at A = 
when ^(A) = 0. 

Then ip is meromorphic on C except possibly at 0. However, because ip is the 
exponent of the opposite of a subordinator, we have (see page 73 of [B96J) 

lim — >} — 00, 01, 

A^0;SR(A)<0 1p(\) 

then 

li m A.-t^-= lim A^(A) = 0. 

A^0;fi(A)<0 tp(\) A^0;Sff(A)<0 

Also, tp is continuous by the right at the point 0, then lim>,_j.o;SR(A)>o ^V'(^) = an d 
we deduce that ip is again holomorphic at 0. 

In conclusion, ip is meromorphic on C and its poles, necessarily in the half plane 
{A; K(A) < 0}, are the same as those of <p (because of the identity ipip = (p and 
because ip is holomorphic on {A; 3?(A) < 0}). Therefore, these poles are the —jj. 
Now, write 
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where \& is an entire function. Notice the limit, valid for any exponent of a subor- 
dinator ( see proposition 2 of Chapter 1 of }B96j ) 



l im ^l e[ o, + oo[ 

A— >+oo A 



Then 



4>M 



(2.1) *(A) = <3(A" +1 ) (3?(A)^+oo) 
And, since 

* w n;(A + 7j ) ,ha) ' 

</>~(A) and ^|^y are 0(A 2 ) on {5R(A) < 0} (see proof of proposition 12. ip and 

is bounded (by ^ * . ) on {5ft(A) < —1}, because it is the Laplace transform of a 

measure supported by ] — oo, 0] (see chapter 3 of |B96) for example) and degP < 
deg Q + 2 , we deduce that 

(2.2) #(A) = 0(A"+ 2 ) (K(A) -> -oo) 

When joining property (|2.1|) and ()2.2j) . ^ is a polynomial with degree at most n+1. 
Moreover, since ^p^- is bounded for A — > +oo, then has a degree at least n. 
Thus 'J is a polynomial and its degree is n or n + 1 . 
If \I/ is a polynomial of degree n+1 then 

lim e]0, +oo[ and lim $(A) e]0, +oo] 

A— >- — oo A A— > — oo 

(this last limit is valid for any exponent of the opposite of a subordinator). Since 
(j) — ipqp^ we obtain 

A— >■— oo A 

And, if "3/ is a polynomial of degree n, then 

lim ^(A)e]0,+co[ an d um e ] _ 00j o], 

A— )■— oo A— — oo A 

(this last limit is valid for any exponent of the opposite of a subordinator) and we 
obtain 

A— > — oo A 

Finaly, ip does not vanish on the left half plane {3?(A) < 0} because it is the exponent 
of the opposite of a subordinator and so, the roots of ip, which are necessary on 
the left half plane {Re(X) < 0}, except possibly (because the exponent of a 
subordinator doesn't vanish in the half plane {i?e(A) > 0}), are the one of <fi. Let 
us study when is a root of if). 

The function tp has a left derivative at (possibly equal to — oo) and the function 
ip has a derivative at because it is a rational function and is not a pole. Then, 
we can write (here, the symbol ' has to be understood as a left derivative ) 



0(0) = ^(o$(o) 
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0'(O) = ip'(0).ip(0) +$(o)v>'(o). 

Since ip(0) e [0,+oo[, tp(0) e [0,+oo[, ^'(0) e [-oo,0[, ip'(0) e]0,+oo] (that 
are obvious consequences of Levy-Kinchin formula for exponent of subordinator or 
opposite of a subordinator), we obtain : 

0(0) = =^> hH°) = and 0'(O) > 0^ or ^(0) = and <f>'(0) < 0^ 

And the proof is finished. □ 

In the sequel we will normalize ip by setting ipoo = 1. 
We will always denote 

— Pi) ~Pli ■ ■ ■ i —fim 
the (possibly equal) roots of ip and 

mi = )J{fc < i;0k = Pi} 

( rrii = for all i's iff the roots — . . . , — f3 m of ip are simple). 

Remark If X is of bounded variations and has a non positive drift then m = n. In 

all other cases, we have m = n + 1, ip is then the exponent of a subordinator with 

positive drift and X "creeps upwards" (see theorem 19 chapter 6 of |B96j ). 

Now, we introduce few notations which will be useful in the sequel, and which are 

related to the negative Wiener- Hopf factor ip. 

First, ip being the exponent of the opposite of a subordinator, the inverse of ip, 4, 
is the Laplace transform of a measure supported by ] — oo,0]. We will denote this 
measure by U(dy), 

-j±-=: [ e-^U{dy) (R(A) < 0) 

WW J]-ao,0] 

Also, for all x g]0,+oo[, we denote by XJy_ x ^(dy) and U]-oo.-x[{dy) the mea- 
sures l[- Xl o](y)U(dy) and 1 ] _ 00 ,- x [(y)U(dy) and by f7[_ a , ](A) and U]-oo,-x[(^) their 
Laplace transform; 



U[- Xl0] (X) := / e- x «U(dy) (A e C) 

J[-x,0] 

#]-»,-*[(*) := / e-^U{dy) (R(A) < 0) 



-oo, — x\ 



Using standard facts about the relation between negative Wiener-Hopf factor ip and 
the fluctuations of the Levy process (see chapter 6 of |B96| ). we obtain the next 
proposition. 

Proposition 2.4. 1) 7/0(0) > or 0'(O) > then then the function |^ is the 

Laplace transform of the measure ip(0)U(dy) which is the distribution of m, the 
minimum of X (m := inf {X t ; < t < 

2) ^(A)C/i_ oc) _ x r(A) is the Laplace transform of the distribution P(Xt x € dy;T x < 
0, (T x :=M{t;X t < -x}). 

Now, using the dual properties settled in this proposition which involves the positive 
Wiener-Hopf factor ip instead of the negative one ip and using the explicit form of 
ip given in theorem 12. 3[ we obtain, (see also Mordecki [LM08 theorem 2.2), 
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Corollary 2.5. 1) If </>(0) > or <j)'(0) < then the distribution of the maximum, 
M :={sup;X t ;ie [0,C[} is 

m m . 

P(M € dy) = aoSo(dy) + V a.^—e-^dy, 

i—l 

where the coefficients a,i are given by the rational expansion : 
Iff =1 (l + A/ 7J ) _ A 

n^(i + A/A) ao + ^(A + /3 l r. +1 

2) The distribution of the "over shoot" X T * — x on T x < ( (T x = mf{t; X t > x}) 
is the following 

n n . 

P(X T « -x G dy;T x < C) = coM^^y) +^c,(x)^e^n y> od2/, 
where the coefficients a,- (a;) are given by the rational expension : 



Cj(a;) 



n^ x (i + A/ft) 



E 



CO 



Corollary [23] and Proposition 12.41 give us a foretaste of the explicit results that we 
obtain in the sequel : Most of the explicit distributions that we will obtain will be of 
the preceding form or of convolution of measures of that form; that is, a combination 
of exponential or gamma distributions with explicit coefficients, possibly restricted 
to an interval, and convoluted with (a possible restriction of) the measure U{dy). 
The next proposition introduces Laplace transforms of distributions that will be 
involved in the sequel. The easy proof is left to the reader. 



Proposition 2.6. 1) 

mi(x, A) := 



r e (^)^ K0 ](-^)-^, iO ,(A) 



-dpy 



A + /3 



is the Laplace transform of the measure 

(±ye{o,x]y m 'e~ PtV dy) * (l ye[ _ Xt0] U(dy)) I l ye [-xfi] 



2) 



n j(x, A) 



Qn 3 



-Xxjj 



] — OO , — X 



[(A) - ei'tJy 



c[ (-7) 



is the Laplace transform of the measure 



3) 



A + 7 

(y n te-™l y>0 )*(U(dy-x)l y<0 > ) I, „ 



ij(X) = 



1 



(-di) n i la + 7 

is the Laplace transform of the measure 

y^e-^n y>0 dy 
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3. Some matrices of coefficients 
Recall first all of the notations already introduced. 

-7i>---)-7n 

are the poles of (/> (and ip) located on the half plane {A, 5ft(A) < 0}, and 

rij = §{k < j;7 fe = 7j} 



— Pi, . . . , —p m 

are the roots of <fi (and ip) located on the half plane {A, 5?(A) < 0} together with 
is 0(0) = and 0'(O) > 0. and 

m = ${k < i;/3 k = Pi} 

We will denote by 1 and c the line and the column : 

1= (l nj =o;l< 3 < n). 



c = (l m4 =o; 1 < i < mf. 
The 111,(3:, A) and nj(x, A) and ij(X) being defined in the previous proposition 
we denote by m(ar, A) the column 

m(x, A) = (v3.i(x, A); 1 < i < m)* 

and by n(x, A) the line 

n(x,X) = (nj(x, A); 1 <j < n) 

and by i(A) the line 

i(A) = (i J (A);l < j <n) 

Also, define the line v(x), 



v(x) := 
and w(x) be the column, 

w(x) := 



(Z^7 e7Xf/ ]-— [(-7) 



1 < i < n , 



7=7j 



(9'' 



U\~ x o}(— P) :l<i<m 
.(-dp)" 1 ' 1 ' JV J/J=ft ~ ~ 

Proposition 3.1. If m = n + 1. i/iera i/ie Zzmii 

fi(-x):= lim (-A)e- A -[7 ] _ 00 ,_ :!;[ (A) 

A— ► — 00 

exists in [0, +oo[ /or a/Z x 6]0, +oo[. 

This proposition will be proved incidently here as a consequence of our computa- 
tions. In fact, it can be shown that the function u is a density of the measure 
l y< oU(dy) and we will see that this density is of bounded variations (see remark 4 
after corollary I5.2j) 

If m = n + 1, we will denote w'(x) the column ^w^(x); 1 < i < to^ , with 

r) mi r i 



w<(x) := 



(-a/3)' 
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Notice that if u is continuous, w^(x) is the left derivative of Wj(x), hence the 
notation. 

Let now W(x) be the matrix (m, n) : 



where, 



W„l.r: : 



W(x) = (W id (x); 1 < i < m, 1 < j < n) 

d mi d n * \e^U [ _ xfi] (-p) + e^U ] _ 00 ^ x[ (- 1 ) 



{-dP) m *{-d~l) n i 



0-1 



When m = n + 1, denote W(x) the square matrix 



W(x) = [ w(x) W(x) ] 
We will show later the next proposition 

Let W be any matrix and let f be any line, the dimension of f being equal to the 
number of columns of W, W[ will be the matrix obtained from W by replacing 
the i-th line by f . And W, will denote the matrix obtained from W by just taking 
it's i-th line off. 

Similarly, let e be a column whose dimension is equal to the number of lines of W, 
W 3 e will refer to the matrix obtained from W by substituting e to it's j-th column. 

4. The distribution of (S^,X^,I^) 

Theorem 4.1. The potential kernel of the triple (S t ,X t ,I t ) is characterized by the 
following identities. For all Ai, A 2 € C, x G]0, +oo[, 



E( [ e- x ^l St - h < x e- x ^-^dt) = E( [ 
Jo Jo 



- A 2 5 t ^ 



s t -i t < x e 



Ut) 



U[-x,o\iM) v(x) 
m(x,Ai) W(x) 



1 



1 



|W(x)| 

U[- x ,o]{\i) u(-x) v(x) 
m(x, Ai) w(x) W(x) 



|W(x)| 
- x * x m{x,\ 2 ) W(x) 



\W(x)\ 

Notations. Let, if m = n, 

lVx,o](0) v(x) 
m(x,0) W(x) 

And, if m = n + 1, 



a(x) := 



a(x) 



\W(x)\ 

1 1 

m(x, 0) W(x) 



if m = n 



if m = n + 1 



r(x) := |W(x) 



a(x) 



U[- x ,o](0) u{-x) v(x) 
m(x, 0) w(x) W(x) 



a(x) := I m(x, 0) W(x) 



r(x) := |W(x)| 



Corollary 4.2. If </>(0) > 0, then the law of the triple (S^,X^,I^) is characterized 
as follow : 



0(0) 



d(x)a(x) 
r 2 (x) 
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Given S{ — < x, the random variable Xq — is independent of Iq and has the 
same distribution as S( and 

1 / m \ 

P(I C G dy|S c -/ c < x) = "a}")' 01 [^oU(dy)-[l ye[ ^ 0] U(dyMj2 h(x)y m 'e-^l ye%x] dy\J 



1 / 

P(S C G ^|5 C -/ C < x) = -'^(ao^(dy)-[l I , e[ o,x]^-»)]*E a *( !,: )y m ' e ~ Al ' 1 ve[o^]dtf] 
with, if m = n, 

da = a = 1 



. , . \Wj {x) (x)\ \W\(x)\ . 

And, if m = n + 1, 

ao = 1 ao = 
_ |» W W->-^)l W | 

|W(x)|| |W(3J)| 

The identities are still true when 0(0) = and when -P^c — — x )i 
P(S C G dy\S c -I(<x) P(J C G dy|5 c - J c < a) 

r+co 

are replaced by / P(S t — It < x)dt, 
Jo 

J +o °P(S t -I t <x,S t £dy))dt J+ 00 P(S t -I t <x,I t edy)dt 

J +o ° P(S t -J t < x)dt / +o ° P(S t -I t < x)dt 

Remark Let </>o be a proper (not killed) Levy exponent and apply the formula of 
theorem to <f> = 4>q + q when q vary from to +oo, let P q be the corresponding 
distribution, that is the distribution of the Levy process X, with Levy exponent 
<f>o, killed at an independent exponential time of rate q. Clearly we have 

rC _ r+°° 



/ e- XlI n St - It < x e- x ^-^dt) = E ( / e-^lst-z^e-^-^e-^dt) 
Jo Jo 

When taking the inverse Laplace transform over q, Ai and A2, we find the measure 

Po(-Tt G dy, Xt — i* G dz, S t -I t < x)dt 

Then the previous theorem and corollary give a characterization of the distribution 
of (St, Xt, It) at any time T. 

As an application, one can compute the price of a double barrier knock-out call 
option under exponential Levy model for the asset Y, Y t = Y e Xt , when the under- 
lying Levy process X satisfy the assumption 12.11 that is, an european call option 
with strike K which is desactivated if the asset Y goes down a barrier value a or 
up a value b before the time of maturity T. It is clear that the price is given by the 
expression : 

E ((loe XT - i«:)l XT>log ^l/ T >iog^-ls T <iog^) 
Then, this price can be deduced from the distribution of the triple (St, Xt, It)- 
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5. Behavior of the process (A, S) at time V x = inf{f, X t - S t < -x} 
Theorem 5.1. 



E Cxp(-/i 2 S'y x -^i(Xy x -5yJ)ly x<C ) = 1p{fXi 



U]-oo, 


-x[(Ml) 


e^ lX .n(x, fii) 


w 


(x) 


W(x) 




1 


1(^2) 






w(x) 


W(x) 





if m = n 



and 



Notations. Denote if m = n, 
1 -i(0) 



HlU]-oo-xl(Hl) -?/]-oo,-x[(Ml) e 


.n(x, ^i) 


w'(x) 


w(x) 


W(x) 




A*2 -1 i(/x 2 ) 






w'(x) w(x) W(x) 





b(x) :- 



w(x) W(x) 

And, if m = n + 1, 

-1 i(0) 
w'(x) w(x) W(x 

Corollary 5.2. 



c(x) := lim V>(A) 



if m = n + 1 



U]-oo,-x[W -n(x,A) 
w(x) W(x) 



&(x) := 



c(x) := lim V>(A) 

A->0 



Af/]-cx>,-x[(A) -Uj-oo ,-x[(A) n(x, A) 
w'(x) w(x) W(x) 



P(K < 



c(x) 
6(x)' 



The random variables Sv x and A14 — Sv x are independent conditionally to V x < Q 
and the distribution of Xy x — Sy x on V x < ( is characterized by the identity, 

U(dz) * P(AV X - S Vx G dz; V4 < C) 

n 

= [6(x)]-^[l 2< _ x t/(^)] * [c.i^^ + coCx)^^)-^^^)^^-^^^^^ 

i 

And, 

m 

P(Sy x e dy; V X <C) = c(x) (Wo + £ h(x)y m ^ e-^l y>0 dy) . 

i=l 

Where, if m = n, 

c_i=0 c (x) = -1 Cj(x) = - W(x)| w(x)] j = 0,...,n 
And if m = n + 1, then 



c_i = 1 c,(x) = W(x)|; +1 



w'(x)] 



i = 0, ...,n 



The fli(x), the exponents rrii(x) := fj{fc < i; (3k(x) = (3i(x)} and the coefficients 
&i(x) are given by the following fractional expansion 



(A + 7j)"i+ lJ 
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-i(A) 
W(a?) 



if m 



if m = n + 1 



A -1 i(A) 

w'(a;) w(a;) W(x) 

Remark 

1) Notice that, when t/>(0) = 0, that means that X does not drift to 
2, chapter 6 of [B96] ) then 

P(X T ^ < C) - P(X Tx < +co) = 1, 

thus 



-oo (see lemma 



B[ (A) = E(e 



-AX 7 



; T x < +oo) ->• 1 A -» 



and we get 



lim ib(X)n.j(x, A) = lim 7 — ^— — 

A->-0 rW JV ' ' A^O {-d~f) n i 



^(A)e- A *l/ ] _ 0Qi _ B[ (A) - V'(A)e^C/ ] _ 00! _, [ (- 7 ) 



A + 7 



7=73 



a « 3 



(_a 7 )« 3 - 



7=7j 



= ii(0) 



Thus, we see that b(x) = c(x) when ^(O) = and, according to the previous 
corollary, we obtain that P(V X < +00) = 1 when X does not drift to +00: this 
means that the heights of the excursions out of the set {t, X t = St} are unbounded 
a.s. This can be obtained by elementary trajectorial arguments. 

2) In order to obtain a direct characterization of the distribution P(A^ — Sy x S 
dz; V x < C) instead of the characterization of the convoluted distribution U(dy) * 
P(Xy :c — SV X £ dz; V x < C)j it is necessary to introduce the drift, the Levy mea- 
sure and the killing rate associated to the exponent tp. We leave this part to the 
interested reader. 

3) Due to the lake of memory of exponential distributions and the fact that 

the distribution of Xy x given Xy m > is, like the distribution of SV X as stated in 
the preceding corollary, a linear combination of the functions y mi ^e~ l3i ^ v l y - > ct- 

4) It is clear from the characterisation of the distribution of Xy x — Sy x when m = 
n + 1 that the convolution l z< - x U(dz) * 5'(dz) is a signed measure, that means 
that the measure U (dy) has a density, and this density, which is the fonction u 
appearing in lemma UTTl is of bounded variations. 



6. Behavior of the process (I, X) at time V x 
Theorem 6.1. 



E(e 



-li 2 (X v x-I v x)\ 





i( M2 ) 






c W(x) 






1 n(x,fii) 




c W(x) 



m£{t; X t -I t > x} 



if 771 







1 










c w(a;) 


W(x) 






1 e-wft 


[(Mi) n 






c 


~w(x) 


W(x) 



if m 
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Notations If m = n 

c(x) -- 

If m = n + 1 , 



i(0) 
c W(x) 



b(x) = lim ip(X) 



1 n(x, A) 
c W(x) 



c{x) 



1 i(0) 
c w(x) W(x) 



b(x) = hmJ(X) 



1 e Ax f7]_oc ,-x[(-A) n(z,A) 
c w(x) W(x) 



One can notice again that b(x) = c(x) when ip(0) = 0. 
Corollary 6.2. 



P(y* < C) 



c(x) 
&(*) 



and /yx and X yx — I V x are independent conditionally to V x < ( 



\Xv—Iv G s+dy|V x < C) = [c( a; )]- 1 .( C o(a ; )5o(dy)+^c j (a ; )^ e -^n !/> ody) 



i=i 

+ 00 

r(x) 

with 



P(Jy. G dz; < C) = ^ [u(dz)] * M r (dz) 



n=0 



fi x (dz) = [l z<Q .U(dz - x)] * [co(x)S (dz) + ^2cj(x)z n * e ^ z l z>0 dz] l z<0 , 
and, if m = n, 



c (x) = and Cj(x) = j^j^ |;<)1 = 1 " 

If m = n + 1, 



c,(z) = ^g-M for ,=0,...,n 
|W(a:)| 

7. Behavior of the process at time U x = mf{t, S t — h> x}. 
Theorem 7.1. 

P(U X < C, X Ux - luj - P(U X < Q,X Ux - SuJ - 

Given the event {U x < Ci^(/ X = Iu x }> &u x and Xjj x — Sjj x are independent. The 
random variable S[/ x has the same distribution as given S^— Iq < x if (f)(0) > 0, 

JTP(S t G dz,S t -I t < x)dt 

and has the distribution — ^ if 6(0) = 0. The random 

f™P(S t -I t <x)dt J Vy ' 

variable Xjj x — Su x has the same distribution as Xy x — Sy x given V x < £. 

Given the event {U x < (;X U:c = Su x }, Iy % and X Ux — I Ux are independent. The 

random variable I Ux has the same distribution as Iq given Sq — Iq < x if 0(0) > 7 

f n °° P(I t £ dy, S t -I t < x)dt 

and has the distribution — — — : if (f>(0) = 0. The random 

J Q P(S t - I t < x)dt 

variable Xjj x — Ijj x has the same distribution as Xy x — Iy given V x < (. 
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8. Behavior of the process at time = mf{t;X t $ [—a, 6]} 
Theorem 8.1. We have, if (f)(0) > 

c(b-y)P(I c edy,S c <b) 



P(X n =s T >X<0= J 



a,o] a ( b - y) 0(°) 



And 



T>(Y T s t\ f c ( a + x ) F (k > -«> £ dx ) 

J[o,b] a(a + x) 0(0) 

k 

0(0) 0(0) 



7/0(0) — 0, we replace in last formulas o^rf Pf-rc^^Scg'fa) &y respec- 

tively 

r + rxi r+oo 

/ P(7 t e dj/,5 t < anrf / P(S t £ dx, -I t > -a)dt 

Jo Jo 

Given {X Tb = S T b,T^ < £}, the distribution of I T b is the distribution of 1^ given 

a (Jo °° P (It edy,S t <b)dt 

J Q + °°P(S t <b,-I t >-a)dt 
0(0) = 0. The distribution of X T b — I T b given I T b = y is the one of X V b+ y CJIV6TI 
yb+y < £ 

Given {X T b — i-p,,T* < (}, the distribution of S^b is the distribution of given 

J + °°P(St £ dx,-I t > -a)dt 

J +O ° P(S* < b; -I t > -a)dt 
0(0) = 0. The distribution of X T b — S T b — x is the one of Xy a+X given 

V a +x < C- 



{S( < b,I{ > —a} if 0(0) > and the distribution — - — if 



{S^ > —cl,Iq < b} if 0(0) > and the distribution +QO — — if 



Remarks 1) We can deduce from this theorem that the distribution of X T b — b given 
X T b > b is a linear combination of distribution of the distributions y nj e~~ fjV l y> ody. 
That could easily be deduced from the lack of memory of the exponential distribu- 
tion. The previous result allows also to compute a closed form of the coefficients. 
2) As in the remark following theorem 14.11 let 0o be a proper (not killed) Levy 
exponent and apply the formula of theorem to = 0o + q when q vary from to 
+oo, let P = P g be the corresponding distribution, that is the distribution of the 
Levy process X, with Levy exponent 0o, killed at an independent exponential time 
of rate q. The probabilities P q (X T b = S T i,T% < () and P q (X T b = I T *,T% < C) 
are q times the Laplace transforms of the functions of t, 

P Q (X T b = S T b , T h a < t) and P [X T t = S T b , T h a < t) . 

On the other hand, the price of a double barrier with rebate, say rebate r + (resp. 
r~) if the option exceed the value b (resp. goes down the value a) before time T. 



(Denote a — — log y-, b = log y-) is given by the next formula 



E((Y e x - - K)l XT>log ^l lT> _~. ST<i )+r+P(X T! = 5 t| ,T| < T) 

+r-P(X T - k =I T i,Ti<T) 

We have seen that the first term can be computed with the help of theorem 14.11 
and now we see that the two other terms can be computed by an inverse Laplace 
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transform in time of P q {X T i = S T i , T| < Q and P g {X T i = I T i ,T~ < () which are 
characterized in the preceding theorem. 

9. Recall of the main result of |F10j and its application with the 

Assumption 12.11 

Proposition 9.1. There exist unique functions A{x, X), A(x, X), B[x, A) and B(x, X), 
C(x,X) and C{x,\) such that A(x, A) and A(x, A), are defined for A £ C and en- 
tire; B(x, A) and C{x, A) are defined, continuous on {5ft(A) > 0} and holomorphic 
on {5ft(A) > 0}; B{x, A) and C(x, A) are defined and continuous on {K(A) < 0}, 
and holomorphic on {5i(A) < 0}, satisfying the following identities, 

P( t e-^l St _ It<x e-^ x t-^dt) = P( / V^l^^e-M*-^) 
Jo Jo 

= i(a;,A 1 )yl(a;,A 2 ) 
E(ex P (-/i 2 Sv, - n(X Vm S Vm )\.v.<t)) = §7^4 

E(exp( -y,il v * - /j, 2 (X v * - Iv°>)lv*<c)) = ^i^'^l 
V / B{x,in) 

%l){X)A{x, A) -)■ 1 B(x, A) - ip(X) for A +oo 

5(x,A)~^(A) for A^-oo 

Proof. The existence of the 6 functions A, A, B, B, C, C has been settled in the 
general setting of a Levy process in proposition 6.1 and proposition 6.2 of |F10] , 
the behavior of the functions is given in theorem 4.2 of |F10J . The uniqueness is 
obvious. 

□ 

More over, we know from theorem 4.2 of jFlO] (or we can deduce them from propo- 
sition inH]), the behavior of the six functions as settled in next proposition 

Proposition 9.2. The functions of X, A(x,X), e~ Xx A(x, A) , e Xx C(x,X) 
are bounded on {Re(X) > 0}. 

The functions A(x,X), e Xx A(x,X), e~ Xx C(x, A) are bounded on {Re(X) < 

0}. 

For Sft(A) — > +oo e Xx C{x, X) i/ ip(X) is bounded, and for K(A) — > —oo, 
e~ Xx C(x, A) -> ififiW * s bounded. 

Remark The property "ip(X) is bounded" means that is the exponent of a 
compound Poisson process and this is equivalent to the property (see chapter 6 of 
|B96j that the first time the Levy process X visits ]0, +oo[ is non zero a.s. (here it 
is also equivalent to the condition m = n). 

In the next lemma, we state that the determination of these six functions reduces 
to the computation of four polynomials. 
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Proposition 9.3. There exist unique polynomials ofX, Pi(x,X), Qi(x,X), P2{x,\), 
Q 2 (x, A), with 

degPi = n degP 2 = m degQi < m — 1 degQ 2 < m — 1, 

such that 

A( x A) = p i(^ A )+Qi( I - A ) e ' A ^N,Q]( A ) 



i(a;,A) = 



B(x,A) 

S(x,A)=Vi(A) 



nr(A + A) 

Q 2 (x,X)e Xx + P 2 (x,X)U^ xM (X) 

nf(A + A) 

P 2 (x,A) 

n?(A + 7,-) 

(Pi (a, A) - Qi(x, A)e- A ^ ] _ 00 ,_ X[ (A)) 



C(x,A) 



nf(A + 7i) 

e- AK gi(a;,A) 



n»(A + 7j ) 

(P 2 (z, A)^ ] _ 00 ,_ a:[ (A) - Q 2 (z, A)e A *) 



n?(A + 7 ,-) 

Proof, of proposition \9.3\ , proposition \3.1[ We have seen in [F10] theorem 4.2 that 
the matrix M(x, A) defined as follows 

(-3$) «w<* 

satisfies the identity 

(9.1) Af+OMu) = M-(x,iu) ( J ) 

Moreover, det M(x, A) = 1, and M is invertible. On the other hand, one can easily 
check that the matrix 

" (a ' A) = ( ^(A) ? ) ( nr(A o +7j) nr(A +ft) ) fOT *W > 



^(A).D' ] _ 0O ,_ x[ (A) f7[-a;,o](A) / V n?(A + 7j ) 

satisfies the same identity 19.11 

Thus, the product N(x, X)M~ 1 (x, A) is entire. One can see easily that, because of 
the boundary conditions fulfilled by the components of M(x, A) given in proposition 
193) that this matrix N(x, X)M~ 1 (x, A) is of the form : 

«(»,A)M-H»,A) . ( ^ lfeA » ) - *(*,A) 

Where Pi, Qi, P 2 and Q 2 are polynomials of A. Then 

det R = dot N = n™(A + /3i)n?(A + 7j -), 
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and we obtain 

M(x,X) = R- 1 (x,X).N{x,X) 

Pi(x,X) e- Xx Q x {x,X) \ ( 1 
e Xx Q 2 (x,X) P 2 (x,X) J X { LV*,o](A) 1 





x i for B(A) > 

V U n«(A+ 73 ) / 

Px(a;,A) e-^Qi^A) \ / $(A) 1 

e^Qa^.A) P 2 (x,A) ^ v -^(A).C/ ] _ 00 ,_ :c[ (A) E/ h *,o](A) 

x | IW+i) ° 1 ] for SR( A ) < 
\ U n-(A+ft) / 
Developping each term of this matrix product, we obtain the identities of the the- 
orem. Now, according to proposition 19.21 e~~ Xx C(x, A) is bounded on {5ft(A) > 0} 
and goes to if tp is bounded (that is if m = n), that gives that degQi < m— 1. 
Futhermore, the equivalence of the proposition 19.11 B(x 7 X) ~ ip(X) gives that 
P 2 (x,X) = X m for A +oo. Since tp(X)A(x, A) 1, we deduce that Fi(A) - A" 
for A — > +oo. 

Also e~ Xx C(x, A) is bounded for K(A) — > +oo and goes to if V>(A) is bounded. 
This is equivalent to say that 

P 2 (x,A)e- A ^ 1 _ oc ^ z[ (A)-g 2 ( a ;,A) 

— — ■ — — ■ — ; > it MX) — > — oo 

H?=i(A + 7j) 

In other words, 

(9.2) P 2 (x, X)e- Xx U ] _ 00 ^ x[ (X) - Q 2 (x, A) = o(A") 

If to = n then P 2 (x, A) ~ A" and we have clearly that e~ ,- x [{X) — > These 

two facts put together with property (|9.2[) give us that Q 2 (x, A) = o(A") and so, 

deg Q 2 (x, A) < n — 1 = to — 1 
If m = n + 1, then P 2 {x, X) ~ A™ +1 and P 2 {x, X)e- Xx U ] _ OCi _ x[ (X) = o(A™ +1 ) then, 
with property (|9.2p . we deduce that Q2(:Ej A) = o(A™ +1 ), in other words, 

deg Q 2 (x, A) < n = to — 1 

The uniqueness of the 4 polynomials come from the uniqueness of the 6 functions. 
This finish the proof of proposition 19.31 

Now, when to = n + 1 and when looking more precisely the preceding identities, 
since degQ 2 (x, A) < n, we have that n ^ 2 ^' A ^ has a finite limit. Since 

P 2 (x, X)e- Xx U ] _ 00t _ x[ (X) - Q 2 (x, A) 

n? =1 (A + 7i ) 

we deduce that 



->0 A^-c 



P 2 (x,X)e- Xx U, 00 _ x[ (X) 

— r has a finite limit when A — > — oo 

n? =1 (A + 7i ) 

Since P 2 (x, A) - A™ +1 and IT/ =1 (A + 7j ) - A", we obtain that 
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Ae Ui—aoi—xnX) has a finite limit (— u(— x)) when A — > — oo. That proves propo- 
sition O Now since p ^ x ^ e ~^-j^iW-Q^ x ^ goe s to 0( for we deduce that 



A™ 



-u(—x). We state this property in the next lemma. 



□ 



Lemma 9.4. 



Q 2 (^,A) 
lim = — u(— x) 

|A|-H-oo A™ V ' 



10. Computation of A and A and proof of theorem 14.11 and corollary 



Proposition 10.1. 

A{x,X) 
A(x,X) = 



1 



1 



e-* x m(x,\) W(x) 



\W(x)\ 

e- Xx m{x, A) W(x) 
|W(x)| 



Proof. Recall the equation of proposition 19.1 

1 



A(x, A) 



when m = n, 



when m = n + 1 



P 1 (x,X)+Q 1 (x,X)e~ Xx U l ^ 0] (X) 



nnA + ft 

Since A(x, A) is an entire function and because degPi = n and degQi < m — 1 (see 
proposition I9.3|) . then ^4(cc, A) is necessarily of the next form 



A(x, A) = a Q (x) - 2J di(a:) 



(-a/3)' 



,3 + A 



do 



{x)-^a t {x)e Xx ini(x,X) 



for some coefhcients ai{x) 1 i = 0, 1, . . . , m 
Thus, 

y m 

Pi (a, A) =n™(A + /3 i )(ao(^) -X>(a!) 



(9'' 



L ( a; ,A)=nr(A + ft) PTa^) — 



(-a/3)' 

gnu 



A + /3 



(-a/3)' 



A + /3 



Now using property Pi(x, A) ~ A™ of lemma [^751 one obtain 



(10.1) 



ao(x) = 1 if m = n 



m 

(10.2) do(x) — and ^""^ (x)wj (x) — — 1 if m = n + 1 

i=i 
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On the other hand, take the identity 



B(x,X) 



P 1 {x, A) - Q x {x, X)e- Ax U ] _ 00) _ x[ (X) 



n(A + 7,0' 

and replace P\ and Q\ by the above expressions, we get 



ni"(A + A) 
n?(A + 7j -) 



d" 



(-apy 



[e^Uj-^j-P) + e- Aa £/ ] _ 00 ,_ a;[ (A) , 
/3 + A 



The function B(x, A) is holomorphic on the left half plan, in particular at points 
A = — 7j, this implies that for every j = 1, . . . , n, 



d 1 ^ ^ . . 9 

i=l 



(5A) 
In other words 



/3 + A 



9" 



/3=/9i,A=- 7j 



(<9A) t 



[oo(x)]a=- 



y^ y ai(x)Wi,j(x) = l nj=0 if 



»=i 



And 



aj(x)Wj j (x) =0 if rn = n + 1 



These n equations added to equation (| 10 . 2[) if m = n + 1, form a linear system of 
to equations with variables a* (a:), i = 1, . . . , m. It is easy to be convinced that each 
solution of this systems brings a new couple of polynomials Pi(x, A) and Q\(x, A). 
The uniquiness of such a couple leads to the uniquiness of this solution. Thus the 
system is a Cramer System and the determinant of the matrix W(x) for to = n 
(resp. of the matrix W(x) for to = n + 1) does not vanish. Finally we obtain, 
For m = n, 



a (x) = 1 Oj(x) 



|W(x) 



for i = 1, 



and 



A(x, A) = a (x) — aj(x)e Ax mi(x,A) = 
l 

For to = n + 1 , 

a (x) = a;(x) • 

and 



1 1 

er Xx va.(x, A) W(x) 



|Wi(x)| 



|W(x)| 

A(x, A) = - V] a i (x)e~ A;r m i (x, A) = 



|W(x)| 

for i = 1, . . . , to, 

e- Aa; m(x,A) W(x) 
|W(x)| 



□ 



Proposition 10.2. 



i(x,A) = 



£7[_ X)0 ](Ai) v(x) 
m(x, Ai) W(x) 



if m = n 
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m(x, Ai) w(x) W(x) 



A(x,\) = 

\W(x)\ 

Proof. Take the equation of proposition 19.11 



if m = n + 1 



A{x, A) = A)e A * + P 2 (z, A)[/ [ _ :C) o ] (A) 

Since A) is an entire function of A and -P 2 (A) ~ A m , degQ2 < m — 1 according 
to proposition 19.31 then A(x,X) is necessarily of the next form 



A(x,X) = U[- x ,o](X) +J2^( X ) 



d m * r E/[_ a>0] (A) - e^ x U hxfl] (-(3) 



(-0/3)' 



X + (3 



= U[- xfl ]{X) - '^2a i (x)m i (x, A) 
l 

For some coefficients cii(x), i = 1, . . . , m. Thus, 



Q 2 (x,X) = -nf(A + Pi)(£2&i(x) 



(-dfsy 



^&[_ X ,0](-A) 



771 

p 2 (^,A)=nf(A + A)(i + X]^( a 



A + /3 
1 



(-dfjy 



X + f3 



13=13, 



Take now the identity of proposition 19.11 

d{x ' A) = nix + jj) ( P2(x ' V e ~ XxiJ ]-oo ,-*[W - Q*( x > A )) 

C(x, A) is holomorphic on {Re(X < 0}. Then for all j = 1, ... , n, 

rP2(x,X)e- Xx U ] ^ 00 ^ x[ (X)-Q 2 (x,X)- 



(dX) n i 



nf(A + ft) 







In this last equation, when replacing P2{x, A) and Q2{x, X) by their above expres- 
sions, we obtain the next system, for j = 1, . . . , n, 



Joxy^ 



e- Xx U } ^ x[ (X) 



. , d n i d r ' 

i=l 



(8X) n J (-dfsy 

In other words , 



- Xx 



A + /3 



/9=/9 < ,A=-7j 



(10.3) 



This a system of n equations where the variables are the m terms di(x). 
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When m = n + 1, we have the property of lemma 19^41 *^ 2 j^' A ^ —¥ —u{—x) when 
A — > —oo then 



^a,(x) 



d 1 ' 



In other words 
(10.4) 



^ a,(x)wj(x) = u(- 



We have already seen in the previous proof, that detW(x) ^ if m = n and 
detW(x) ^ if to = n + 1. Then the system (| 10 . 3[) augmented by the equation 
(|10.4[) if to = n + 1, is a Cramer system; and the solution is given as follows. If 

to = n, 



Oi(x) = 



\w [ : (x)] (x)\ 

\W{x)\ 



and 



A(x, A) = U^ X)0 ](\) - ^2 Oi(x)mi(x, A) = 
i 

If m = n + 1 , 



C/[-x,o](A) v(x) 
m(x, A) W(x) 



|W(x)| 



|W(x) 



&i(x) 

|W(x) 

m 

I(x,A) = ^_ x ,o](A) - ^a l (x)e Aa; m i (x,A) = 



[u(-a;),v(x)]| 



U[- X , ](X) u(-x) v(x) 
e Aa: m(a;. A) w(x) W(x) 



|W(x)| 



□ 



Proof of theorem 14.11 and corollary 14.21 

Theorem 14.11 follows from ropositions 19.11 110.11 and 110.21 

After that, corollary 14.21 becomes obvious by the Laplace inversion of the functions 
rrii(x, A) given in proposition 12.61 

11. Computation of C(x, A) and Z?(x, A) and proof of theorems 15.11 and 

COROLLARY 15.21 



Proposition 11.1. 



C7(x,A) =V>(A) 



U]-oo.- x[ ( x ) -e Xx n(x,X) 
w(x) W(x) 



|W(x) 



B(x,X) 



1 -i(A) 
w(x) W(x) 



|W(x) 



if m = n 



C(x,A) = ^(A) 



XU]-oo^x[(X) -L>]_oo,-^(A) e Xx n(x,X) 
w'(x) w(x) W(x) 

|W(x)| 
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B(x,X) = 



A -1 i(A) 
w'(x) w(x) W(x) 



if m = n + 1 



|W(aO| 

Proof. We take the expression of C(x, A) of proposition 19.3 

C(x, A) = n n ( ^ ) (-0a(«. A ) eA * + A)«7 Hoo ,_ x[ (A) J 
with dcg P2 = m and deg Q2 < m — 1. 

Because C(a:, A) is holomorphic on {5R(A) < 0}, the function e~ Xx C(x, A) is in the 
next form, 

e- Ax C(x,A) =^(A).(c_ 1 (a;)Ae- A ^ ] _^_ K[ (A) + c (x) e - A ^ ] _ 00 ,^ [ (A)+do(a:) 



3 = 1 

In other words, 



e' Ax ^]-oo,-.[(A) - e^^_ 0O| _ a[ (- 7 ) 



A + 7 



7=7j 



e- Ax C(x,A) 



n 

^(A)(c_i(x)A.e- AK ^ ] _ M ,_ K[ (A)+c (x)e- AK L> ] _ oo ,_ K[ (A)+d (^)-^Cj(x)n J (x,A) 



for some coefficients Cj(x), j = — 1, 0, 1, . . . ,n and a coefficient do (a?) 
Thus 



P 2 (x,X) = n(A + 7 J )fc_i(a:)A + co(x) -^Cj(z) 



a ™ 3 



1 



A + 7 



7=7j 



And 



71 

Q 2 (x,X) =n(A + 7 J -)(-*(x)-£ej(a 



0"' re^t/ ] „ 00 ,_ x[ (- 7 ) 



A + 7 



7=7i 



(-fry)*- 

Take the equation of proposition 19.11 

e- x *A(x,\) = — L — (Q 2 (x, A) + ft (g, X)e- Xx U [ - xfi] (A)) . 

The function A(x, A) is an entire function, then we obtain the equation for every 
i = 1, . . . , m, 



0" 



(dxy 



2 (x,X) + P 2 (x,X)e-^U hxM (X) 



n?(A + 7i) 







Taking the expressions of P2 and Q2 upwards, we obtain the next equations, 



(dX) m 

n 
3=1 



A2_i(a;)e Xx U^ xfi] (X) + c (x)e Xx U { _ xfi] (X) - d Q (x) 



(<9A) m * (-fry) B * 



A + 7 



A=-ft,7=7j 
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In other words, 
(11.1) 



co(a;)Wj(a;) + ^ Cj(x)W i)j (x) 



(d\y 



We can compute the coefficients C-i(x), co(x) and do(x) by using the assertions of 
lemma l§751 More precisely, if m = n, we have P2(x, A) ~ A" and degQ2 < n— lj 
thus we obtain 



(11.2) 



c-i(x) = d (x) = cq(x) = 1 



If m = n + 1 then A) ~ A m and according to proposition 

[— u(— x)]X n , thus we obtain 



Q2(a;, A) 



(11.3) do(x) = u(—x) c_i(x) = 1 

Then, if m = n, using the system (jll.ll) and equation (|11.2p we get the new system 

n 

Cj- (x) Wjj (#) = -Wi(i) for i = 1, . . . ,n 

Thus 



Cj(x) 



|W? , Jx 
|W(s)| 



And 



C(x,A)=^(A).(c/ ] _ 00 ,_ :c[ (A)-e Aa .^c i (x)n i (a : ,A)) 

i=i 

^]-oo,-x[(A) -e^nfoA) 
w(a;) W(x) 



V'(A) 



|W(ar)| 



And 



P 2 (x,A) 

n(A + 7j ) 



71 Oyj ■ -i Tl 

= 1 -E^p^-^U, = 1 - 



1 -i(A) 
w(z) W(x) 



|W(z)| 

Similarely if m = n + 1, we obtain from system (|ll.ll) and equation (|11.3I) . the new 
Cramer system 



Cp (x)Wj (x) + ^ Cj (x)W it j (x) = 



d" 



3 = 1 



(-9/3)' 



^ x U { _ xfi] {-p) + u{-x)\^ =v?l(x) 



Then (for j = 0, 1, . . . , n), 
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Cj(x) 





x)\ 




W(x) 





And 



C(x, A) = ^(A).(A.0' ] _ oo ,_ B[ (A)+fi(-a:) 


e AK +c (a;)t)' ] _ OOi _ :!; 


n 

j'=i 




A^i-oo ,-x[(A) +u(-a;)e A:r -Oi^ x[ (A) 


e Aa n(x, A) 




w'(x) 


w(x) 


W(x) 






W(x) 



B(x, A) = T ^ 2 X (X ' *\ = A + go (g) - ^ c^x 



n(A + 7j ) uw V (-*r) ni 

n 

= A + c (x) - ^Cj^ijiX) 



A + 7 



7=7i 





j=i 




A 


-1 


i(A) 


w'(x) 


w(x) 


W(x) 




W(x) 





□ 



Proof of theorem 15.11 and corollary 15.21 

The proof of theorem 15.11 follows from the above proposition and proposition 19.1 
Corollary 15.21 follow then with the help the inverse Laplace transform of the func- 
tions nAx, A) given in proposition 12.61 

12. Computation of the functions C(x,X) and B(x, A) and proof of 

THEOREM 16. II AND COROLLARY 16. 21 



Proposition 12.1. 



C(x, A) =e 



- A.l' 



i(A) 
c W(x) 



|W(x)| 

C(x,X) = e' Xx 



1 n(x, A) 
c W(x) 



|W(x)| 



if m = n 



1 i(A) 
c w(x) W(.t) 



B(x,\)=4>(\) 



|W(x)| 

e - A ^ ] - 00 ,_ ;c[ (A) n(x,A) 
c w(x) W(x) 



|W(x) 



i/ m = n + 1 



Proof. According to proposition 19.31 we have the identity 

^(A) 



B(x,A) = 



.(Pi(x,X) - Qi(x, X)e- Xx U ] ^ oc ^ xl (\)) 



with 



n?(A+ 7j -) 

deg Pi (A) - A™ and deg Q x < m - 1 
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Since B(x, A) is holomorphic on {K(A) < 0}, we can write B(x, A) in the following 
form, 

0»* re^U^.^i-^-e-^U^^X)- 



B(x,X)=i>(X)i l+b (x)e- Xx U ] ^ oo ^ x[ (X)+Y,b 3 ^)- 



j=l 



-dlY 



A + 7 



In other words, 

n 

B(x, A) = $(A). (l + b (x)e- Xx U ] _ oo> _ x[ (X) - £ 6, (x)n 3 (x, A)) 

3=1 

for some coefficients bj{x), j = — 1, 0, 1, . . . , n and a coefficient do(x). 
Thus the polynomial Pi and Q\ are given by the next expressions 



P 1 (x,A)=nT(A + 7i) (l + £>(x) 



(9" 



3=1 



(-d 7 )"i 



A + 7 



7=7j 



A) = HT(A + 7i ) -6o(x) + £ &,(*) 



<9" 



3=1 



(_3 7 )« 3 la + 7 



1 



7=7j 



On the other hand, according to proposition IHHJ we have 



A(x, A) = W n ( ^ +(3i) fa*, A) + Qi(x, X)e- Xx U [ _ xfi] (X)) 

Since A(x, A) is an entire function of A, then it is holomorphic at A = —ft; Thus 
we have the equations for i = 1, . . . , m. 

d m * rA(^A)+g 1 ( a :,A)e-^[/ [ _^ 0] (A)- 



that is, 



n?(A + 7i ) 



A=-/3< 



0. 



6 (a;) 



e Aa: ^[-x,o](A) 



A=-/3« 



3=1 

In other words, 



(dX) m *(-d-/) n i 



(dxy 

e 7X C/]-oo,-,[(-7) + e- Aa; t/ [ - :c ,o](A) 



A + 7 



7=73i^=- Pi 



— l mi =0 



7=7j 



(12.1) h (a;)w i (a;) + ^ ^ (a;) W i;i (a;) = l m , =0 

3=1 

Moreover, if to = n, since degQi(x, A) < n — 1, we have 



(12.2) 



So (aO = 



We have already seen that the determinants |W(x)| ( if m = n) and |W(x)| (if 
to = n + 1), then the system (|12.1[) augmented with (|12.2I) if to = n, form a Cramer 
system whose solutions are given by the next expressions, 
If m — n, 



|W c (a;) 

b (x) = bj{x) = 



\W(x)\ 



for j = l,...,n 
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Thus 



1 n(x,X) 
c W(x) 



\W(x)\ 



and the expression of C(x, A) given in proposition [973] gives 



If m = n + 1, we have for j = 0, . . . , n, 

bj(x) = 



(E 



Xx 



i(A) 
c W(ai) 



f A + 7j 



|W(x)| 



And 



A) = V-(A) (l - SoCasJe-^Oi-oc-xtCA) - J2 ( x ' A ) 



V-(A) 



1 e- Xx Uy 


oo-x[(A) n(x,A) 


c w 


(x) 


W(x) 




W(x) 





and then, 



CO, A) 



Qi(x,A)e- Aa: 

' n?(A + 7j ) 



f A + 7j 



-A.r 



1 i(A) 
c w(a:) W(a;) 



W(z) 



□ 



Proof of corollary 16.21 This corollary becomes obvious when using proposi- 
tions 19.11 and the previous one and when using the Laplace thransform inversion of 
nj(x,X) and ij(A) given in proposition 12.61 



13. Proof of theorems 17. II and 18. II 

It has been proven in |F10j (propositions 6.3. and 6.4) that the distributions of the 
triple (X, I, S) at time U x and at time T\ are also characterized by the 6 functions 
A, A, B, B, C, C. Then, we translate these two propositions with reference to the 
explicit forms of the functions obtained previously. 
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